Version: 11 March 2009, 11:00 EDT 



The Chandra COSMOS Survey, I: 
Overview and Point Source Catalog 

Martin Elvis^, Francesca Civano^, Cristian Vignali^'^, Simonetta Puccetti^'^, Fabrizio 
Fiore^, Nico Cappelluti^, T. L. Aldcroft^, Antonella Fruscione^, G. Zamorani^, Andrea 
Comastri^, Marcella Brusa^''', Roberto Gilli^, Takamitsu Miyaji^'^, Francesco Damiani^", 

Anton Koekemoer^^, Alexis Finoguenov^'^, Hermann Brunner^, CM. Urry^^, John 
Silverman^^, Vincenzo Mainieri^^, Guenther Hasinger^'^^, Richard Griffiths^^, Marcella 
CaroUo^^, Heng Hao^, Luigi Guzzo^^, Andrew Blain^*^, Daniela Calzetti^^, C. Carilli^", 
Peter Capak^^, Stefano Ettori^, Giuseppina Fabbiano^, Chris Impey^^, Simon Lilly^^, 
Bahram Mobasher^^, Michael Rich^^, Mara Salvato^*^, D.B. Sanders^^, Eva Schinnerer^^, N. 
Scoville -^^y Patrick Shopbell^^, James E. Taylor^'', Yoshiaki Taniguchi^^, Marta Volonteri^^ 



- 2 - 



^Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 USA 
■^INAF— Osscrvatorio Astronomico di Bologna, Via Ranzani 1, 1-40127 Bologna, Italy 
^Dipartimcnto di Astronomia, Univcrsita dcgli Studi di Bologna, Via Ranzani 1, 1-40127 Bologna, Italy 
^ASI Science Data Center, via Galileo Galilei, 00044 Frascati Italy 

^INAF— Osservatorio astronomico di Roma, Via Frascati 33, 00040 Monteporzio Catonc, Italy 

^Max-Planck-Institute fiir Extraterrestrische Physik, Postfach 1312, 85741, Garching bei Miinchen, Ger- 
many 

^University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 

^Instituto de Astronomia, Universidad Nacional Autonoma de Mexico, Ensenada, Mexico (mailing ad- 
dress: PO Box 439027, San Ysidro, CA, 92143-9027, USA) 

^Center for Astrophysics and Space Sciences, University of California San Diego, Code 0424, 9500 Gilman 
Drive, La JoUa, CA 92093, USA 

^•^INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 1-90134 Palermo, Italy 

^^Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218 USA 

^^Department of Physics and Yale Center for Astronomy & Astrophysics, Yale University, P.O. Box 208121, 
New Haven, CT 06520-8121, USA 

^^Department of Physics, Eidgenssische Technische Hochschule-Zurich, CH-8093 Zurich, Switzerland 

"ESQ, Karl-Schwarschild-Strasse 2, D-85748 Garching, Germany 

^^Max-Planck-Institute fiir Plasmaphysik, Boltzmannstrasse 2, D-85748 Garching bei Miinchen 

^^Departmcnt of Physics, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA 15213 USA 
^^INAF— Osservatorio Astronomico di Brcra, via Bianchi 46, 23807 Moratc, Italy 

i^CaUfornia Institute of Technology, MC 105-24, 1200 East Cahfornia Boulevard, Pasadena, CA 91125 UA 
^^Department of Astronomy, University of Massachusetts, Amherst, MA 01003 USA 
2°National Radio Astronomy Observatory, PO Box O, Socorro NM 87801 USA 
^^Spitzer Science Center, 314-6 Caltech, Pasadena, CA 91125 USA 

^^Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721 USA 
^^Department of Physics and Astronomy, University of California, Riverside, CA 92521 USA 
^^Department of Physics and Astronomy, University of California, Los Angeles, CA 90095 USA 
^^Institute for Astronomy, University of Hawaii, 2680 Woodlawn Dr., Honolulu, HI 96822 USA 
^^Max-Planck-Institut fiir Astronomie, Konigstuhl 17, Heidelberg D-69117, Germany 
^''Department of Physics and Astronomy, University of Waterloo, Waterloo, Ontario, N2L 3G1, Canada 



-3- 



ABSTRACT 

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra 
program that has imaged the central 0.5 sq.deg of the COSMOS field (centered 
at 10'*, +02°) with an effective exposure of ~160 ksec, and an outer 0.4 sq.deg. 
area with an effective exposure of ~80 ksec. The limiting source detection depths 
are LgxlO"^*" erg cm"^ s"^ in the Soft (0.5-2 keV) band, 7.3x10"^^ erg cm"^ s~^ 
in the Hard (2-10 keV) band, and 5.7x 10"^^ erg cm'^ s"^ in the FuU (0.5-10 keV) 
band. Here we describe the strategy, design and execution of the C-COSMOS 
survey, and present the catalog of 1761 point sources detected at a probability 
of being spurious of <2xl0-^ (1655 in the Full, 1340 in the Soft, and 1017 in 
the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing 
positions with the ACIS-I imager, a remarkably uniform (±12%) exposure across 
the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower fiux 
limit. The widely different PSFs obtained in each exposure at each point in 
the field required a novel source detection method, because of the overlapping 
tiling strategy, which is described in a companion paper. This method produced 
rehable sources down to a 7-12 counts, as verified by the resulting logN-logS 
curve, with sub-arcsecond positions, enabling optical and infrared identifications 
of virtually all sources, as reported in a second companion paper. The full catalog 
is described here in detail and is available on-line. 

Subject headings: surveys - catalogs - X-rays:general - cosmology:observations - 
(galaxies:) quasars: general - galaxies: evolution 

1. Introduction 

The co-evolution of galaxies and quasars or active galactic nuclei (AGNs) has been vig- 
orously pursued both observationally and theoretically for a decade, ever since the discovery 
that the mass of the central black hole is tightly correlated both with the luminosity (Magor- 
rian et al. 1998; Marconi & Hunt 2003) and the velocity dispersion of the spheroid {Mbh-o' 
relation; Ferrarese & Merrit 2000; Gebhardt et al. 2000; Tremaine et al. 2002). Tackling 
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this large subject requires the study of both galaxies and AGNs back to the epoch when 
both were growing rapidly, i.e. z~l-3, requiring deep observations across many wavelengths, 
from radio through the infrared, optical and ultraviolet, to the X-rays. At the same time, 
the wide range of cosmic density and the rapid changes in this large scale structure (LSS) 
require wide field observations that sample the Universe at close to their true fractions. 

The Cosmic Evolution Survey (COSMOS, Scoville et al. 2007a) is a deep and wide 
extragalactic survey designed to have sufficient area to overcome most cosmic variance, which 
requires sampling regions some 50 Mpc on a side (Fig.l; Scoville et al. 2007a), and with 
sufficient depth to sample the z = 1 - 3 galaxy and AGN population. The contiguous 2 sq.deg 
COSMOS field samples a volume of -6x10^ Mpc^ at z = 0.5 - 1 (Wright 2006). This is 
~10% of the volume imaged by the Sloan Digital Sky Survey (SDSS) in the local (z<0.1) 
universe (5.7xlO'^Mpc^, 8000 sq.deg, DR5^). COSMOS is a region of low, uniform. Galactic 
obscuration (E(B-V)~0.02 mag, Nh (2.7xl02°cm-^ Dickey & Lockman 1990). COSMOS 
is likely to be the largest survey of this type for the next decade. 

The location of the COSMOS area near the equator (10'^, +02°) allows all major and 
future facilities^ (notably EVLA, ALMA, and SKA) to target this region down to faint 
hmits (Scoville et al. 2007a). Space-based imaging has been undertaken in the F840W 

(~ i-band) with Hubble Space Telescope (HST, Scoville et al. 2007b), in the 3.5 fim-70fim 
infrared using Spitzer IRAC and MIPS (Sanders et al. 2007), in the UV using GALEX 
(Zamojski et al. 2007), and in 0.5-10 keV X-rays with XMM (Hasinger et al. 2007, Cappclluti 
et al. 2007). Ground-based imaging spans the radio (1400 MHz VLA, Schinncrcr et al. 
2007), the near-lR with CTIO and KPNO (Capak et al. 2007) and CFHT (McCracken et 
al. 2009, in preparation), the optical to AB~26-27 with Subaru in 21 bands (Taniguchi 
et al. 2007). Finally, large dedicated ground-based spectroscopy programs in the optical 
with Magellan/IMACS (Trump et al. 2007), and VLT/VIMOS (Lilly et al. 2008) are well 
underway. 

This wealth of data has resulted in an initial 15-band photometric catalog of ~10^ ob- 
jects (Capak et al. 2007) from which photometric redshifts good to <3% for z<1.2 and r <24 
have been derived (Mobasher et al. 2007). Recently, more photometric bands have been 
added, resulting in improved photo-z's for the galaxy population accurate to Az/(l+z)<l% 
(Ilbert et al. 2009) and to Az/(l+z)~2% for the AGN population (Salvato et al. 2009). 

We have undertaken the Chandra-COSMOS survey (C-COSMOS) to cover the central 
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Fig. 1. — Left: Map of the COSMOS field sliowing the coverage at various wavelengths: the IRAC 
3.6fj,m mosaic is the background image; Blue solid= C/iandra, blue dashed= C/iandra deep; black 
polygon = HST; red solid = Subaru, CFHT, zCOSMOS bright; red dashed = zCOSMOS deep; 
black dashed = XMM and VLA. The Spitzer MIPS observations cover an area 2 times larger. 
Right: Large scale structure seen in galaxy distributions in the COSMOS field (Scoville et al. 
2007c), ranging in size from 1-20 Mpc, based on photo-z's accurate to Blue is centered on 

z=0.35, Green on z=0.55, and Red on z=0.75, each with Az=0.05. The C-COSMOS field outline 
is shown as the white tilted square, with the dashed line delineating the high exposure area as in 
the left panel. A scale showing 10 Mpc at the three redshifts is shown at the top. In both panels 
North is up, East is to the left. 

0.9 sq.deg region of the COSMOS field (Fig.l, left), containing a wide range of cosmic 
overdensity (Fig.l, right), with the ACIS-I CCD imager (Garmire et al. 2003) on board the 
Chandra X-ray Observatory (Weisskopf et al. 2002). The survey took 1.8 Msec of Chandra 
observing time (~21 days) and was the largest guest observer program approved in a single 
AO at the time it was undertaken (2006 November - 2007 June). C-COSMOS employed a 
series of 36 heavily overlapped ACIS-I 50 ksec pointings to give an exposure of ~160 ksec over 
the inner area to a depth of ~1.9xl0^^^erg cm~^s^^ (0.5-2 keV), providing an unprecedented 
combination of contiguous area and depth in the X-ray band. This overlapping tiling strategy 
gives highly uniform exposure, and so a well-defined flux limit. 

Several of the deepest COSMOS surveys are now concentrating on this same central 
sub-field of COSMOS: the z-COSMOS Deep spectroscopic survey (to B~25, Lilly et al. 
2007), the deep VLA survey (6 /xJy rms, Schinnerer et al. 2009, in preparation), and several 
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millimeter and sub-millimeter surveys (MAMBO, Bertoldi et al. 2007 and AzTEC, Scott 
et al. 2008). GALEX has observed the central field deeply (Zamojksi ct al. 2007) and is 
currently monitoring this area. The Ultra- VISTA survey will undertake a deep yJHK survey 
of the central 1.5 sq.deg, half of which will be surveyed to the unprecedented limits of 
--26 AB mag (Arnaboldi et al. 2007). 

By going for large area rather than extreme depth, most of the C-COSMOS sources are 
sufficiently bright to be detected in the rest of the pan-chromatic COSMOS data set, allow- 
ing rapid identifications (Civano et al. 2009) and determination of their multi-wavelength 
properties (e.g. Elvis et al. 2009, in preparation). On the other hand, C-COSMOS is suf- 
ficiently deep that significant numbers of normal and starburst galaxies with luminosity of 
10^^ erg s~^ can be detected up to z ~0.9, a redshift depth comparable with that of the 
galaxy redshift surveys in the COSMOS field (Taniguchi et al. 2007, Lilly et al. 2007). 
Adding the Chandra coverage to the COSMOS survey adds a valuable resource for the study 
of the co-evolution of black holes and their host galaxies, of the SEDs of faint quasars and 
active galactic nuclei, and the evolution of galaxies. 

The summed image of the entire C-COSMOS field is shown in Fig. 2 where colors have 
been mapped to X-ray bands. 

This is the first of three papers presenting the basic results of the C-COSMOS survey 
over the whole field. Paper I (this paper) reports on the strategy, design and execution of 
the C-COSMOS survey, and present the catalog of 1761 point-fike X-ray sources detected in 
C-COSMOS; Paper II (Puccetti et al. 2009) presents the details of the simulations carried 

out to optimize the source detection method; Paper III (Civano ct al. 2009) presents the 
identification of the X-ray sources with optical and infrared counterparts. We conclude by 
listing the primary science objectives foreseen for the C-COSMOS data. Papers on several 
of these topics are in preparation. 

We assume a ACDM cosmology with Ho=70 km s"\ 0^=0.27, O^ac=0.73. 

2. The Chandra COSMOS Strategy 

For C-COSMOS we have developed a strategy that uses ~50% overlapping tiling of the 
16.9x16.9 arcmin ACIS-I fields. This tifing produces a remarkably uniform sensitivity in the 
central part of the field, and a well-defined fiux limit with a sharp cut-off (Fig.3; for details 
on the generation of sensitivity maps see §7 in Paper II). This approach also ensures that 
the area with HPD<2" is maximized, so that the unique Chandra high resolution imaging 
(van Speybroeck et al. 2002) can be exploited fully, albeit with 1/4 of the exposure time. 



Fig. 2. — X-ray image of tlie Chandra COSMOS field, optimized to show point sources with a 
wide variety of X-ray colors. A scale showing 30 arcmin (the approximate diameter of the full 
moon) is shown for comparison. North is at the top; East is to the left. The full angular resolution 
of Chandra is not well represented in this image as, in order to display the point sources clearly, 
the original image has been smoothed with a sharp gaussian with radius equal to 2.9", and added 
to an image of the field smoothed with a wide gaussian with radius equal to 4.4". X-ray 'colors' 
are mapped so that red is the 0.5-2 keV band, green is the 2-4.5 keV band. Blue is the 4.5-7 keV 
band, and each energy band was smoothed in the same way. Selected prominent clusters have been 
adaptively smoothed for display (red extended shapes). 
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Table 1: C-COSMOS flux limits and corresponding XMM-COSMOS flux limits. 



Band 


C-COSMOS(lim)« 


C-COSMOS(IogN-logS)" 


XMM-COSMOS" 


Soft (0.5-2 kcV) 


1.9 


2.5 


5 


Hard (2-10 kcV) 


7.3 


16 


25 


FuU (0.5-10 keV) 


5.7 







" flux limits are reported in units of 10 ^^erg cm ^s ^, for bands up to 10 keV, but were 



measured only up to 7 keV. (See text for details.) 

The good Chandra point spread function (PSF) resolves sources 2" apart over ~0.7 sq.deg, 
corresponding to 8-16 kpc separations for 2;=0.3-0.9, and locates point sources to <4 kpc 
at any redshift. Thus close mergers can be resolved, and nuclear sources distinguished from 
off-nuclear sources in galaxies (Ultra-luminous X-ray Sources, ULXs, Fabbiano 2006, Lehmer 
et al. 2006, Mainieri et al. 2009, in preparation). 

Point source detection sensitivities were estimated for three standard Chandra bands: 
Soft (S, 0.5-2 keV), Hard (H, 2-10 keV) and Full (F, 0.5-10 keV). Due to the high background 
in the 7-10 keV energy range^, channels above 7 keV were not used for source detection. (See 
§4.2.2 and Paper II for details). The C-COSMOS flux limits in 3 bands are reported in 
Table 1, together with the XMM-COSMOS hmits for comparison: C-COSMOS sensitivity 
is three times below the corresponding flux limits for the XMM-COSMOS survey (dashed 
line; Cappelluti et al. 2009), making them complementary surveys. 

The achieved sensitivity-area curve^(Fig. 3) has a sharp cut-off at low fluxes. 

The C-COSMOS Soft band flux limit corresponds to luminosities of (0.8, 4, 11) x lO'^^erg s^ 
at 2;= (0.3, 0.6, 0.9) respectively, while the Hard band flux limit corresponds to four times 
higher luminosities. Both luminous elliptical galaxies and starbursts often exceed these lu- 
minosities, and starburst galaxies are known to become common (Hornschemeier et al. 2003) 
at these X-ray fluxes. 

The low ACIS background enables stacking analysis, in which counts at the positions 
of known classes of objects, e.g. subsets of the thousands of galaxies with redshifts, are 
co-added to increase the effective exposure time (Brusa et al. 2002; Hornschemeier et al. 
2002; Brandt et al. 2001; Nandra et al. 2002; Fiore et al. 2008a,b). 



^URL: http://cxc.harvard.edU/contrib/niaxini/bg/index.htnil#spec 

^This curve is remarkably close to the predictions from the proposal, reflecting the high accuracy with 
which the requested tiling was executed. 
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Fig. 3.— Area-flux curve for C-COSMOS (red solid line, 0.5-2 keV). The coverage of ECDFS 
(Lehmer et al. 2005; dashed hne), AEGIS-X (Laird et al. 2008; dash-dotted hne), CDFN (Alexan- 
der et al. 2003; magenta short-long dashed line), CDFS (Luo et al. 2008; cyan dotted line) and 
XMM-COSMOS (Cappelluti et al. 2009; black dashed line) are shown for comparison. 
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2.1. Design 

The C-COSMOS tiling scheme (Fig. 4, left panel) covers the central area of the COS- 
MOS field in the most efficient manner that we could devise. A 6x6 raster array of 36 ACIS-I 
pointings (one ACIS pointing field of view is outlined in black in Fig. 4, left), each of 50 ksec 
nominal exposure, were chosen. The center of the array (Table 2) is shghtly offset from the 
center of the COSMOS field to match the z-COSMOS deep field (Lilly et al. 2007). 

The value of the 8.0' offset between pointing centers was chosen to be slightly less than 
the 8.3' size of an ACIS chip (Garmire et al., 2003, Chandra Proposers' Observatory Guide, 
aka POG, 2007^) Table 6.1), so that chip gaps are not co-added to create small scale dips in 
the effective exposure time. 

The inner part of the field was covered by four exposures, to give a total nominal 
exposure of 200 ksec (effective exposure ~160 ksec) over a 42'x42' area (0.5 sq.deg,green 
area in Fig. 4). The outer region has been covered by two observations (blue area) and 
the four corners covered by 1 observation (purple area) . The corners of the outer and inner 
regions are reported in Table 2 clockwise from the top left. 

Sources at a ffux of ~2.0xl0~^^erg cm~^s~^ (0.5-2 keV) have a total of 5-10 summed 
counts in the four exposures, ensuring a good detection, given the low Chandra/ ACIS back- 
ground of ~2 counts/200 ksec over a 2 arcsec radius circle (see §4.2.1). 

The heavily overlapped tiling scheme produces a smooth exposure map that is ffat to 
12% in the central region (see Figure 4, right panel and §4.2.2). 



2.2. Comparison with Other Legacy Surveys 

Chandra observing time has been dedicated to several large legacy surveys: CDF-S 
(Giacconi et al. 2002; Luo et al. 2008), CDF-N (Alexander et al. 2003), ECDF-S (Lehmer 
et al. 2005), AEGIS-X (Nandra et al. 2005; Laird et al. 2008), XBootes (Murray et al. 
2005). These surveys have different emphases in area and depth, so we summarize the 
special features of C-COSMOS here. 

Like all contiguous area surveys, C-COSMOS has significant advantages over non- 
contiguous surveys (e.g. SEXSI, Harrison et al. 2003, Eckart et al. 2006; ChaMP, Kim et al. 



^Chandra X-ray Center publication TD 403.00.010 
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Table 2: Coordinates of the C-COSMOS field, center and corners of the outer and inner 
regions, clockwise from the NE (top left). 



RA 




Dec 




Center 








10'' 00" 


24-^ 


+02° 10' 


55" 


Outer region 






IQh 02™ 


45* 


+02° 26' 


47" 


QQh 59m 


11" 


+02° 46' 


45" 


QQh 57m 


54** 


+01° 53' 


00" 


IQh 01™ 


23* 


+01° 33' 


59" 


Inner region 






^Qh 02™ 


05* 


+02° 21' 


13" 


QQh 59m 


30* 


+02° 35' 


47" 


09'' 58™ 


35* 


+01° 59' 


19" 


IQh Qim 


11* 


+01° 44' 


37" 



2007) , because of the difficulty of getting deep multi- wavelength coverage of non-contiguous 
fields. 

C-COSMOS is neither the deepest (CDFN and CDFS) nor the widest (XBootes) legacy 
Chandra survey. A comparable sensitivity has been reached in the somewhat smaller AEGIS 
field (dot-dashed line in Fig. 3; Laird et al. 2008). C-COSMOS differs from the other surveys 
by having the largest area at fiuxes 0.3 — 1 x 10~^^erg cm~^s~^, and a sharper low flux limit 
cut-off in the area surveyed than most other recent X-ray surveys. The single fleld CDF-S 
and CDF-N have notably shallower roU-offs in their sensitivity curves (magenta and blue 
lines in Fig. 3). 

Hence, to compare the area and depth of C-COSMOS with comparable contiguous 
Chandra and XMM-Newton surveys in a consistent fashion requires a slightly revised mea- 
sure of area and depth. We have used the Area-Flux plot from each survey to derive the 
flux at the point where each survey reaches 80% of the maximum survey area. We plot 
these values in Figure 5 (filled circles) for the Chandra contiguous area surveys (CDFN, 
Alexander et al. 2003; CDFS, Luo et al. 2008; ECDFS, Lehmer et al. 2005; AEGIS- X, 
Laird et al. 2008; XBootes, Murray et al. 2005; ELAIS-N, Manners et al. 2003), and for the 
XMM-Newton contiguous surveys that fill regions of the fiux-area plane (ELAIS-Sl, Puc- 
cetti et al. 2006; XMM-COSMOS, Cappelluti et al. 2009; Lockman-Hole, Brunner et al. 

2008) . The C-COSMOS fiux at 80% of the area covered (0.72 sq.deg) in the Soft band is 
6xl0~^^erg cm~^s~^ 

Compared with other plots of this kind (e.g. Brandt & Hasinger 2005) survey points in 
Figure 5 tend to be moved diagonally toward smaller area and high flux hmits. This shift 
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Fig. 4 — Left: The 'as designed' C-COSMOS tiling for the 36 50 ksec pointings. The thick black 
box (top left) represents one ACIS-I pointing, the thin boxes all the pointings. Different colors 
show areas with different number of overlapping pointings: green - 4 overlapping pointings; blue - 
2 overlapping pointings; purple - 1 pointing. The black bars show roughly the relative dimensions 
of one pointing (~16'), of the inner area with larger exposure (~42'), and of the total field (~56'). 
Raster point (see Table 3 1-1 lies at the top right (NE) and 1-6 lies at the top-left (NW). Right: 
The 'as executed' exposure map for the C-COSMOS survey in the Soft band. The color bar gives 
the achieved effective exposure in units of seconds. 

can be quite large for surveys with shallow slopes at low fluxes in their area-flux limit curves 
(as for examples the deep flelds). This is because the normally quoted area is the maximum 
area of the survey, while the normally quoted flux limit is that of the faintest source in the 
survey, which can be detected only in a much smaller area. 

Curves of constant numbers of sources (for the Soft band) are shown in Fig. 5 following 
the predictions of Gilli et al. (2007) XRB model^. The larger numbers in XBootes and 
the two COSMOS surveys are notable. Some 1000 sources are predicted for C-COSMOS 
above the '80% area' flux limit in the Soft band based on the logN-logS relation of Gilli et al. 



^The curves have been computed using the tool "POrtable Multi Purpose Application for XRB and AGN 
counts" available at the web site http://www.bo.astro.it/~gilli/counts.html. 
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(2007), while 1023 are actually detected. This can be compared with the CDF fields which 
have ~200 sources each. 

A single number does not convey the complete picture, of course. We also show in Fig. 5 
the Area-Flux curve of each survey down to 20% of the area. These curves better explain the 
differences between the surveys, notably between the two CDF deep fields, that are due to 
the different observation's strategy [changing only the roll angle (CDF-N) or also moving the 
centroid (CDF-S)]. The more sensitive, smaller area, parts of each survey add more sources 
than indicated by the dashed black lines, especially for the curves that are closer to vertical. 
For example, the AEGIS-X survey (Laird et al. 2009) has 1032 soft sources, about double 
the number predicted at the 80% point. C-COSMOS, with a flatter flux-area curve, has a 
total of 1340 S band sources, ~30% higher than the 80% area number. 

Each of these surveys has extended mult i- wavelength coverage, but C-COSMOS is the 
only deep and wide X-ray survey fleld selected for both existing deep multi-wavelength 
coverage, and for future legacy value, due to the equatorial location of the COSMOS field. 
The AEGIS field (5=+52°), the CDF-N field (5=+62°) and the XBootes field (5=+35°) are 
all too northerly to be accessible by ALMA or the VLT. The COSMOS field was also selected 
to have low IR cirrus emission, and a lack of bright stars. X-ray or radio sources in the field 
to maximize multi- wavelength coverage. 

C-COSMOS and XMM-COSMOS complement one another by providing large samples 
of sources over a wide fiux range (Fig. 8), while sharing the same extensive multi- wavelength 
data set. XMM-COSMOS provides a larger sample of extended sources, while C-COSMOS 
provides a larger sample of starburst and normal galaxies. 

3. Observations 

A summary of the Chandra ACIS-I C-COSMOS observations as carried out is given in 
Table 3. Primarily because of thermal constraints on spacecraft components (POG, §3.3.3), 
many of the 36 C-COSMOS pointings were scheduled as two or more separate ObsIDs, 
giving 49 C-COSMOS observations in all. The indices X-Y (1-1 through 6-6) describe the 
fleld numbers, where X is an index in RA and Y an index in Dec, with 1-1 being in the top 
right (NE) corner of Fig. 4 (left panel), and 1-6 being in the top left (NW) corner. 

The observations took place in two main blocks: 2006 December - 2007 Jan and 2007 
April- June (Table 3). The flelds were observed at nominal roll angles of 250/70 deg, where 
the visibility of the COSMOS fleld is at maximum (~70%) and the pitch angle is such that 
the constraints are either unrestricted or restricted only to avoid overheating of the charged 
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Fig. 5. — Area-Flux curves for Chandra (red) and XMM-Newton (blue) contiguous X-ray surveys. 
Each survey has been plotted using each sensitivity curve starting from the flux corresponding to 
the area that is 80% of the maximum area for that survey (large points at the top of each curve), to 
the flux corresponding to the 20% of the total area (bottom of each curve). Data were taken from 
the following Chandra references: XBootes - Murray et al. 2005; CDFN - Alexander et al. 2003; 
CDFS - Luo et al. 2008; ECDFS - Lehmer et al. 2005; AEGIS-X - Laird et al. 2008; ELAIS-N 
- Manners et al. 2003; and XMM-Newton references: XMM-COSMOS - Cappelluti et al. 2009; 
Lockman - Brunner et al. 2008; ELAIS-S - Puccetti et al. 2006. The black dashed curves show the 
total number of 0.5-2 keV sources expected based on the logN-logS relation predicted by Gilli et 
al. (2007) at the 80% area point. 
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particle detector (EPHIN). As an equatorial field, the roll angle of the COSMOS field is 

quite stable (55.2-69.6, 248.4-256.2) for periods of ~100 days. As a result, the Chandra 
X-ray Center (CXC) Mission Planning team were able to maintain a tight roll angle range 
of ±6° around the nominal values (Table 3), leading to a highly uniform exposure of the 
whole field. 

The mean effective exposure time per field (not per ObsID) is 46.3 ksec, when only the 
Good Time Intervals (GTls), cleaned of the few high background times (§4.2.1) are used''. 
The maximum exposure is 48.3 ksec and the minimum exposure (excluding a single 37.6 ksec 
exposure for field 2-5, Table 3) is 44.1 ksec. So, with this exception, the range of exposures 
over the fully covered inner region varies by just ±2.0 ksec (4%). 



4. Data Processing 

The data from the 49 obsids were uniformly processed in two phases using the ClAO 
3.4 software tools^ (Fruscione et al. 2006) , the yaxoiP tool and custom versions of the XMM 
SAS detect tool EMLdetect^'^. Standard Level- 1 and Level-2 processing pipeline^^ (AS- 
CDS version 7.6.9) from the CXC were used. In the first processing phase we determined 
astrometric corrections (see below) for each ObsID. These corrections were then applied in 
the second phase where we reprocessed all event data starting with Level- 1 products. 

Data processing involved the following series of steps, as summarized below: 

1. Astrometric corrections (<1.1") to the standard COSMOS frame starting with the 
CXC supplied standard data products (§6); 

2. Baseline data product creation by re-processing all ObsIDs to a standard frame of 
reference using the new astrometry and standard CXC pipelines (§4.2); 



^This is ~93% of the requested exposure, well within the 90% tolerance limit defined for Chandra schedul- 
ing. 

^URLhttp: / / cxc.harvard.edu/ ciao/ 
^http: / / cxc.harvard.edu/contrib / yaxx/ 
^°http://xmni.esac.esa.int/sas/8.0.0/EMLdetect 

^^Pipeline processing levels are explained at URL: http://cxc.harvard.edu/ciao/data/sdp.htnil 
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Table 3: Chandra- COSMOS observation summary 



Field 


Obs. Start 


Exp. Time* 


RA 


Dec 


Roll 


Obs. ID 


(UT) 


(ks) 


(J2000.0) 


(J2000.0) 


(deg) 



7995 


1-1 


2007 


Jun 01, 


03:41 


44.6 


10 


02 


02 


05 


-1-02 


22 


36.46 


248 


4 


7996 


1-2 


2006 


Dec 28, 


11:28 


44.7 


10 


01 


31 


99 


+02 


25 


20.48 


63 


4 


7997 


1-3 


2006 


Dec 30, 


21:10 


44.5 


10 


01 


01 


92 


+02 


28 


04.50 


62 


8 


8494 


1-4 


2006 


Dec 16, 


13:21 


20.2 


10 


00 


31 


85 


+02 


30 


48.52 


66 


4 


8122 


1-4 


2007 


Jan 20, 


10:15 


28.0 


10 


00 


31 


85 


+02 


30 


48.52 


55 


2 


8493 


1-5 


2006 


Dec 12, 


18:07 


19.3 


10 


00 


01 


79 


+02 


33 


32.55 


66 


4 


7998 


1-5 


2007 


Jan 10, 


21:41 


26.9 


10 


00 


01 


79 


+02 


33 


32.55 


63 


2 


8478 


1-6 


2006 


Nov 24, 


10:17 


17.6 


09 


59 


31 


72 


+02 


36 


16.58 


69 


6 


7999 


1-6 


2006 


Nov 25, 


09:24 


29.0 


09 


59 


31 


72 


+02 


36 


16.58 


69 


6 


8000 


2-1 


2007 


May 26 


20:23 


45.2 


10 


01 


51 


10 


+02 


15 


05.52 


253 


2 


8001 


2-2 


2007 


Apr 02, 


03:42 


47.3 


10 


01 


21 


03 


+02 


17 


49.54 


256 


2 


8123 


2-3 


2007 


Apr 07, 


13:40 


48.3 


10 


00 


50 


97 


+02 


20 


33.55 


255 


2 


8002 


2-4 


2006 


Dec 19, 


04:57 


28.5 


10 


00 


20 


90 


+02 


23 


17.58 


65 





8496 


2-4 


2006 


Dec 23, 


12:05 


17.8 


10 


00 


20 


90 


+02 


23 


17.58 


65 





8003 


2-5 


2007 


Apr 02, 


17:53 


37.6 


09 


59 


50 


83 


+02 


26 


01.61 


255 


2 


8004 


2-6 


2006 


Nov 27, 


02:25 


15.3 


09 


59 


20 


76 


+02 


28 


45.64 


68 


6 


8482 


2-6 


2006 


Dec 02, 


09:05 


10.2 


09 


59 


20 


76 


+02 


28 


45.64 


68 


6 


8483 


2-6 


2006 


Dec 04, 


03:02 


21.3 


09 


59 


20 


76 


+02 


28 


45.64 


68 


6 


8005 


3-1 


2007 


Apr 25, 


02:42 


30.8 


10 


01 


40 


15 


+02 


07 


34.57 


255 


2 


8552 


3-1 


2007 


Apr 26, 


09:33 


14.4 


10 


01 


40 


15 


+02 


07 


34.57 


255 


2 


8124 


3-2 


2007 


Apr 08. 


03:42 


31.1 


10 


01 


10 


08 


+02 


10 


18.59 


255 


2 


8549 


3-2 


2007 


May 05 


17:17 


17.2 


10 


01 


10 


08 


+02 


10 


18.59 


255 


2 


8503 


3-3 


2006 


Dec 31, 


10:18 


20.0 


10 


00 


40 


02 


+02 


13 


02.61 


62 


2 


8006 


3-3 


2007 


Jan 01, 


11:48 


25.8 


10 


00 


40 


02 


+02 


13 


02.61 


62 


2 


8007 


3-4 


2006 


Dec 19, 


22:18 


21.1 


10 


00 


09 


95 


+02 


15 


46.64 


64 


2 


8497 


3-4 


2006 


Dec 25, 


01:50 


27.1 


10 


00 


09 


95 


+02 


15 


46.64 


64 


2 


8008 


3-5 


2007 


Jan 02, 


04:39 


45.0 


09 


59 


39 


88 


+02 


18 


30.67 


61 


9 


8009 


3-6 


2007 


Jan 02, 


18:06 


44.8 


09 


59 


09 


81 


+02 


21 


14.70 


61 


8 


8010 


4-1 


2007 


Apr 27, 


18:45 


32.9 


10 


01 


29 


19 


+02 


00 


03.29 


255 


2 


8553 


4-1 


2007 


Apr 29, 


01:02 


14.4 


10 


01 


29 


19 


+02 


00 


03.29 


255 


2 


8011 


4-2 


2007 


Apr 04, 


04:08 


45.8 


10 


00 


59 


13 


+02 


02 


47.30 


255 


2 


8012 


4-3 


2007 


Jan 04, 


05:30 


48.0 


10 


00 


29 


06 


+02 


05 


31.33 


61 


3 


8013 


4-4 


2007 


Jan 04, 


19:44 


46.9 


09 


59 


58 


99 


+02 


08 


15.36 


61 


1 


8014 


4-5 


2007 


Jan 05, 


09:29 


44.2 


09 


59 


28 


92 


+02 


10 


59.38 


60 


9 


8015 




2007 


Jan 07, 


09:53 


44.1 


09 


58 


58 


85 


+02 


13 


43.42 


60 


2 


8550 


5-1 


2007 


Apr 18, 


19:11 


22.7 


10 


01 


18 


25 


+01 


52 


32.34 


255 


2 


8016 


5-1 


2007 


Apr 19, 


20:24 


23.3 


10 


01 


18 


25 


+01 


52 


32.34 


255 


2 


8017 


5-2 


2007 


Apr 04, 


17:55 


45.3 


10 


00 


48 


18 


+01 


55 


16.35 


255 


2 


8018 


5-3 


2007 


Apr 05, 


07:17 


45.8 


10 


00 


18 


11 


+01 


58 


00.38 


255 


2 


8019 


5-4 


2007 


Apr 06, 


23:25 


48.0 


09 


59 


48 


04 


+02 


00 


44.41 


255 


2 


8020 


5-5 


2007 


Apr 09, 


06:12 


47.8 


09 


59 


17 


97 


+02 


03 


28.44 


255 


2 


8021 


5-6 


2007 


Apr 09, 


20:24 


47.3 


09 


58 


47 


90 


+02 


06 


12.48 


255 


2 


8022 


6-1 


2007 


May 10 


23:28 


30.9 


10 


01 


07 


30 


+01 


45 


01.39 


251 


4 


8555 


6-1 


2007 


May 12 


16:06 


16.2 


10 


01 


07 


30 


+01 


45 


01.39 


251 


4 


8023 


6-2 


2007 


Apr 10, 


12:49 


48.3 


10 


00 


37 


24 


+01 


47 


45.41 


255 


2 


8024 


6-3 


2007 


Apr 11, 


21:40 


47.9 


10 


00 


07 


17 


+01 


50 


29.44 


255 


2 


8025 


6-4 


2007 


Apr 12, 


11:57 


47.9 


09 


59 


37 


10 


+01 


53 


13.47 


255 


2 


8026 


6-5 


2007 


Apr 13, 


07:31 


45.8 


09 


59 


07 


03 


+01 


55 


57.49 


255 


2 


8027 


6-6 


2007 


Apr 14, 


13:54 


48.3 


09 


58 


36 


96 


+01 


58 


41.53 


255 


2 



"After GTI and high-background filtering for two affected obsids. Intervals of 8.50 ksec and 2.45 ksec (re- 
spectively) were eliminated from the two affected ObsIDs (8003, 8014). 
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3. Background reduction using high background time filtering (which affects only two 
ObsIDs) (§4.2.1); 

4. Exposure map creation in the three energy bands F, S, and H, using the standard 
CIAO tool sequence (§4.2.2); 

5. Calculation of the sky coverage (i.e. the area covered to a given flux threshold) in the 
three energy bands, F, S and H; 

6. Candidate source detection using a wavelet technique (PWDetect, Damiani et al. 
1997)12; 

7. Selection of reliable sources, with a probability of being spurious <2xl0~^ in at least 
one band, using maximum likelihood fitting (EMLdetect) applied simultaneously to 
each ObsID at the positions of all candidate sources; Puccetti et al. (Paper II) shows 
that EMLdetect reconstructs the input count rate of simulations well, while both 
PWDetect and detector underestimate the input count rate by about 15%; 

8. Reliability checks for all sources using simulations, searches for outliers and visual 
checks (rejected candidate sources were all in the wings of bright source PSFs); 

9. Aperture photometry of reliable sources. At high fluxes the systematic error in the 
PSF, which is intrinsic to the EMLdetect method, becomes larger than the statistical 
error; this systematic error is not present for aperture photometry. 

10. Derivation of reliability and completeness criteria for the source catalog, leading to a 
logN-logS curve that provides an end-to-end check of the source extraction by compar- 
ing with other surveys in the same flux range (§9). 

Steps 1 to 4 are discussed more fully in the following subsections. Complete details of 
the steps from 5 onwards, including details of the simulations and tests, are given in Paper II. 

4.1. Astrometry corrections 

In the first phase we determined accurate astrometric offsets for each ObslD. The good 
absolute astrometry produced by Chandra (0.6"at 90% confidence, POG, §5) is still of the 



^^We compared PWDetect with the CIAO tool wavdetect used by most Chandra deep surveys on a 
subset of C-COSMOS fields, and found no substantive difference in the results; PWDetect is a much faster 
algorithm, due to better memory buffering. 
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order of one ACTS pixel. To avoid a loss of sensitivity, correcting the astrometry to much 
less than one pixel error is needed before merging event files, or stacking. 

To this end, we first produced a list of bright X-ray sources for each of the 49 ObsIDs, 
using the standard CIAO CELLDETECT tool. Starting with the standard ACTS Level-2 
data products, we generated a broad-band exposure map for each ACTS CCD using the 
CIAO^^ tools ASPHIST, MKINSTMAP, and MKEXPMAP. These exposure maps and event files 
were then used as input to a Chandra- adscptcd version of the XMM-SAS tool EMLdetect 
(see next section), with an input source candidate catalog obtained by running the sliding 
cell detection tool eboxdetect with a high threshold. All sources detected with likelihood 
parameter C >10 were compared with the CFHT MegaCam 1-band catalog of the COSMOS 
field (Capak et al. 2007), selecting only the point- like sources with I magnitudes in the range 
18-23. Using this restricted magnitude range minimizes systematic effects introduced by 
bright stars (saturation) and faint background objects (misidentification) , and is appropriate 
for sources in this flux range (Brandt & Hasinger 2005). An optical-X-ray position correlation 
was computed using the likelihood algorithm included in the SAS task EPOSCORR (Cappelluti 
et al. 2007, 2009). This task uses all the possible counterparts of an X-ray source in the 
field to determine the most likely coordinate displacement. This method is independent of 
the actual spectroscopic identifications, but post facto all the identifications have proved to 
be correct (Paper 111). No statistically significant offset in roll was required for any ObsID, 
so the change in roll was set to exactly zero. The systematic offsets between the X-ray and 
the optical positions were always smaller than 1.1 arcsec, with an average shift of ARA= 
0.04" and Adec=0.25". 



4.2. Baseline Data Products 

The second phase of processing brought the 49 Chandra ObsIDs to a common reference 
frame using the offsets derived above, and generated the baseline data products that were 
then used as the starting point in all subsequent C-COSMOS analysis. 

This processing was based on the CIAO thread for creating a new Level-2 event file from 
Level- 1 products^^. First, a new aspect solution for each ObsID was generated to remove the 
astrometric offset for each ObsID derived in above section, using the REPROJECT_aspect 
tool. Then a new bad pixel file was created using ACIS_run_hotpix (see 'background reduc- 



' URL : http : / /cxc . harvard . edu/ciao/ 

http://cxc. harvard. edu/ciao3.4/threads/createL2/ 



ARA (X-I) arcsec 



Fig. 6. — The X-ray to I-band offsets (ARA, Adec) in arcsec for X-ray sources witli a secure 
identification (Civano et al. 2009) after the aspect correction described in the text has been apphed 
(§4.1). The circles encompass 68% (0.56"), 90% (0.81") and 95% (1.41") of the sources with optical 
counterparts and secure identification. Red dots mark sources with less than 50 counts in the Full 
band. 

tion' below). Finally, a new ACIS Level-2 event file was then created for each ObsID using 
the AClS_PROCESS_EVENTS tool, with: (a) the standard ASCA grade set (grades [0, 2, 3, 4, 
and 6], FOG §6.14), (b) pixel randomization turned off, (c) FHA randomization turned on, 
(d) Very-Faint mode processing enabled, and (e) the new aspect solution applied. 

The astrometric corrections were checked using X-ray sources with point-like optical 
counterparts (Civano et al. 2009, Paper III) that were not used to derive the offsets for the 
individual ObsIDs. The residual systematic shift (X-ray - Optical position) is on average 
Aa= -0.1" and Adec=0.08", and the 1 a dispersion is 0.56" (i.e. the radius within which 
68% of sources lie; Fig. 6). We find that 90% of the X-ray positions agree with the identified 
optical/IR counterpart positions to within 1.1". The residual systematic shift is small enough 
that it will not affect the identification of any individual source and is smaller than the 
average X-ray positional error, and therefore has not been used to correct the astrometry 
any further. The good quality of the data provides positions with sub-arcsecond accuracy at 
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off-axis angle <6', in agreement with other Chandra surveys (0.23-1.90" in the CDFS, Luo 
et al. 2008; 0.3-1.67" in AEGIS, Laird et al. 2008). 

4-2.1. Background reduction and cosmic ray afterglow detection 

Intervals of high background were determined by creating a background light curve for 
the ACIS-I CCD events with point sources found by wavdetect in the phase 1 processing 
removed. Only two obsids showed intervals with a significant (> 5-a) deviation from the 
quiescent background level (see Table 3). 

Particular care was taken in the rejection of cosmic-ray afterglows'^. When a cos- 
mic ray hits a CCD pixel a residual charge can remain localized for tens of seconds and 
produce "afterglow events" , that appear to be X-ray events, at one location for several con- 
secutive CCD frame readouts (POG §6.9). To reject cosmic ray afterglows wc used the 
CIAO tool ACIS_RUN_HOTPIx'^ and enabled Very-Faint (VP) mode background processing 
in ACIS_PROCESS_EVENTS. This process was successful as none of the C-COSMOS sources 
subsequently detected have the time localization characteristic of a spurious afterglow source. 
This procedure also gave a 25-30% background reduction in the 0.5-7 keV band. 

The residual background is very stable over the full field of view at ~ 1.8 x 10"''' 
counts/s/pixel or ~2 counts/200 ksec over a 2 arcsec radius circle, which represents the 
typical size of our detection cell across the field. Following Alexander et al. (2003), in which 
the transition between a photon limited and a background limited regime is defined as >3.3 
background counts per detection cell for S/N=3, we conclude that C-COSMOS is photon 
limited for point source detection. 

4-2.2. Exposure Maps and Sensitivity Curve 

We constructed exposure maps using the standard CIAO tool sequence of asphist, 
MKINSTMAP, and MKEXPMAP, for each ObsID on a per-CCD basis, in each of three energy 
bands, S, H, F. 

Figure 4 (right panel) shows a composite image of the effective exposure time (sec) in 
the Soft band. We clearly see the central region with four overlapping pointings, the side 



URL: http:/ /cxc. harvard. odu/ciao/why/aftcrglows. html 
http: / / cxc.harvard.edu/ ciao3.4/ahelp/acisj:un_hotpix.html 
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Fig. 7. — Histogram of the exposure times in the summed C-COSMOS field. The narrow peaks 
he at the 1, 2 and 4 exposure vahics. The broader bases correspond to overlaps caused by slight 
variations in the roll angles of the ObsIDs. 

strips with two observations, and the corners covered by just one pointing. The uniformity 
of the exposure in the central region is shown by the histogram of the exposure times shown 
in Fig. 7. This histogram shows narrow peaks at the 1, 2 and 4 exposure values, which have 
gaussian sigmas of 12.9, 13.6 and 19.3 ksec, respectively, i.e. a 12% spread on the central 
region exposure. The total effective exposure in the inner, 4 exposure, region is ~160 ksec 
at the peak, and ~170 ksec at the mean, in the same region (see Fig. 7). 

The C-COSMOS sky coverage (i.e. the area covered as a function of hmiting sensitivity) 
was computed in the three standard energy bands, F, S, H using the exposure maps, the 
background maps and assuming a spectrum with F = 1.4 and Nh =N^(Galactic). The sky 
coverage in the Soft band is shown in Fig. 3. More details on the Full band and Hard band 
are given in Paper II (§7). 

The main uncertainty in the estimated sky coverage comes from the range of conversion 
factors from count rates to fluxes induced by the variety of intrinsic X-ray spectra in the 
X-ray population, in both power-law slope and intrinsic absorption, at a minimum. More 
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Table 4: Conversion factors from count rates in the Soft, Full and Hard bands (0.5-2, 0.5-7, 
2-7 keV) to fluxes in the same bands for different spectral assumptions", computed with the 
Chandra Cycle 8 response matrices. 



r 




factor^ 




_2 

cm 




Soft Band 




1.4 


Galactic 


i.O 1 


1.7 


Galactic 


1 81 


2.0 


Galactic 


i. / 


1.4 


1022 


9 1 9 


2.0 


1022 


2.15 


Full Band 




1.4 


Galactic 


0.75 


1.7 


Galactic 


0.89 


2.0 


Galactic 


1.04 


1.4 


1022 


0.51 


2.0 


1022 


0.71 


Hard Band 




1.4 


Galactic 


0.38 


1.7 


Galactic 


0.43 


2.0 


Galactic 


0.47 


1.4 


1022 


0.36 


2.0 


1022 


0.45 



" r=1.4, N^f=Galactic used for catalog fluxes, 
conversion factor CF where Flux = Brate/{CP * 10^^)> in units of cts erg"^ cm?. 

complex spectra are surely present. An additional complication is that the average spectral 
properties are a function of the observed flux (Brandt & Hasinger 2005). To estimate this 
uncertainty, we calculated the sky coverage for power law spectra with F =1.4 and 2.0 with 
Galactic N^j, and for absorbed power law spectra with F —lA and 2.0 and N^j = 10^^ cm~^. 

The range of conversion factors, given by PIMMS^^ is a factor 2.0 in the F band, 1.3 in the 
H band and 1.2 in the S band (Table 4). As expected from the large width of the Full band, 
the uncertainty for the Full band is larger than for the Soft and Hard bands. 



^"^URL: http://http://cxc. harvard.edu/toolkit/pimms.jsp 
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5. Point Source Catalog 
5.1. Overview 

In this catalog we report the 1761 sources detected down to a defined threshold in 
at least one band. The threshold was chosen to balance completeness (the fraction of true 
sources detected) against reliability (the fraction of false sources detected) . Paper II describes 
simulations that allowed us to choose a threshold which has a known completeness and 
reliability. We chose a probability threshold of P=2x 10~^, giving 99.8% reliability for 
sources with more than 12 counts and 99.7% reliability for sources with 7 counts. This 
implies ~3-5 spurious F band detections in the full field with > 12 counts and 5 spurious 
detections with > 7. At this threshold, the simulations then show that C-COSMOS is 87.5% 
complete for 12 count sources and 68% complete for 7 count sources. The C-COSMOS false 
source rates are consistent with those of other surveys (e.g. AEGIS-X, Laird et al. 2008) 
once the higher C-COSMOS threshold and larger average source extraction region are taken 
into account (see Paper II, §6, 8,9). 

The Maximum Likehhood statistic detml = —ln{P) = 10.8 for P=2xl0~^, and this 
threshold detml was applied in EMLdetect. The numbers of source detections at or 
above detml^lQ.% are fisted in the left column of Table 5. Cross-matching the sources with 
detml>lQ.S in the three bands gives a total of 1761 sources. There are numerous sources with 
detml>10.8 in fewer than three bands. In these cases we can search for significant fiux in 
the other bands to a 100 times higher P, as the area being searched is now 100 times smaller 
than the whole survey area (for a 5" cross- match radius). This corresponds to a threshold 
detml =6. In the right hand column of Table 5 we give the numbers of sources detected 
in each band having 6< detml <10.8. Table 6 reports the numbers of catalog sources at or 
above detml — 10.8 in 3 bands, 2 bands, or in only one band. (In this last case the sources 
must have detml>10.S in order to have been selected at all.) 

Almost a thousand (946) XMM-COSMOS sources have also been observed by Chandra 
with an exposure larger than 30 ksec (Cappelluti et al. 2009), and 876 are present in the 
C-COSMOS catalog. Only 70 sources are not present in the Chandra catalog, while 24 
XMM-COSMOS sources have been resolved into two separate sources (Brusa et al. 2009; 
Paper III) due to the better Chandra PSF. Of the 70 sources not recovered by Chandra, 
more than half are in regions with low exposure (between 30 and 50 ksec) as, for example, in 
small gaps of low exposure (Fig. 4). The remainder are either sources with only hard XMM 
detections or, after a visual inspection, they are found to be spurious XMM sources, in 
agreement with the expected fraction of spurious sources. C-COSMOS and XMM-COSMOS 
combine to give a total of ~2800 unique COSMOS X-ray sources. The distribution of X-ray 
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Table 5: Number of sources detected in each band at the two adopted thresholds. 



Band 


detml >10.8 


6< detml <10.8 


Full (F) 


1655 


71 


Soft (S) 


1340 


88 


Hard (H) 


1017 


165 



Table 6: Number of sources with detml >10.8 in at least one band. 



Bands 


Number of sources 


F+S+H 


922 


F+S 


474 


F+H 


257 


F 


73 


S 


32 


H 


3 


Total 


1761 



fluxes for the C-COSMOS sources in the Soft and Hard bands is shown in Figure 8. For 

comparison, we also show the flux distribution of CDFN (dotted line), CDFS (dot-dashed 
line) and XMM-COSMOS detected sources (dashed line). The Chandra and XMM-Newton 
surveys are complementary in that, together, they span almost 3 orders of magnitude in 
X-ray flux, and have over 100 soft band (and over 50 hard band) sources per 0.16 dex bin 
over about 1.5 orders of magnitude in flux. The well-defined cut-off in source numbers 
at faint fluxes, which reflects the tight exposure time distribution (Fig. 7), is significantly 
different from the relatively flat distribution of CDFN (dotted line) and CDFS (dot-dashed 
fine) source fluxes (Fig.8). 

The complete catalog contains source positions and source count rates, exposure times, 
signal-to-noise ratio (S/N), and fluxes in the three bands and hardness ratios (see next 
section). The catalog is ordered with the sources detected in the Full band flrst, followed by 
those detected in the Soft band only and by those detected in the Hard band only. 

The resulting catalog is available in the ApJ on-hne version and on the ' Chandra COS- 
MOS Survey' website^^. Supporting data products (including images, event flies and 



' littp://chandracosmos.cfa. harvard. edu/reports/analysis/20090310_TA_source_catalog_2.1/ 
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Fig. 8. — The distribution of fluxes for sources detected in the Soft and Hard band (continuous 
line) compared with the flux distribution of CDFN sources (dotted line), CDFS (dot-dashed line) 
and XMM-COSMOS sources (dashed line). Sources with upper limit have not been included in 
this figure. 

exposure maps) will be available at the 'Chandra COSMOS Survey' website and at IRSA^^. 
At the Chandra COSMOS Survey it will also be possible to browse a database that includes 
'postage stamps' of the X-ray data for each source, along with the multiwavelength optical 
and infrared data, including the I-band, K-band and Spitzer 3.6/zm (Band 1) images used in 
Paper III to identify the sources. 

5.2. Catalog Description 

The EMLdetect procedure was run on the three bands: Soft, Hard and Full. In order 
to be consistent with other results in literature, count rates estimated in the 2-7 keV and 
0.5-7 keV energy bands were extrapolated into 2-10 keV and 0.5-10 keV fluxes, respectively. 




URL: http://irsa.ipac.caltech.edu/data/COSMOS/ 
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using a spectral slope, r=1.4. We also report the number counts obtained from aperture 
photometry (see Paper II). 

Table 7 gives the columns of the catalog of the 1761 X-ray sources. A more detailed 
description of each column is reported below: 

- Column 1: Chandra source name, following the standard lAU convention with the 
prefix "CXOC" for 'Chandra X-ray Observatory COSMOS' survey 

- Column 2: Source number. Sources are hsted in order of detection: first those detected 
in the Full band with detml >10.8, followed by those detected in the soft band only 
and by those detected in the Hard band only. 

- Column 3-4- Right Ascension and Declination in the J2000 coordinate system. 

- Column 5: Positional error + "^Dec) computed using the following equation 
PoSerror = P S Fradiu.-i / ^/S whcrc S Is thc mimbcr of net source counts, after the sub- 
traction of thc background, in a circular region of radius corresponding to thc 50% 
encircled energy in the field where the source is at the lowest off-axis angle (Paper II). 

- Column 6-7: Count rate and count rate error in the Full band (0.5-7 keV). These are 
effective count rates that would apply if thc source had been observed at the aim point 
in every pointing. I.e. computed by dividing the best fit counts for each source by 
the effective exposure time at the position of each source (the effective exposure time 
includes corrections for vignetting, dither, bad pixels and spatially-dependent quantum 
efficiency). The count rate error at 68% confidence level was computed using the 

equation error == , where where Cg are the source counts estimated 

by EMLdetect, corrected to an area including 90% of the PSF^°, B are the background 
counts evaluated from the background rate (counts/pixel) estimated by EMLdetect 
multiplied for an area of radius R^, which is the mean of the radii, correspondig to 
90% enclosed counts fraction (ECF) of each observation, weighted by the observation 
exposure relative to the total exposure, and T is the vignetting corrected exposure 
time at the position of the source from the exposure maps. We use a=0.5, to allow 
for uncertainties in the background, which is computed through the EMLDETECT 
procedure (see Paper II for more details). 

- Column 8-9: Full band 0.5-10 keV fiuxes and errors were computed converting count 
rates to fiuxes using the following formula: Flux — Brate/{CF * 10^^), where Bj-ate is 



' http : / /cxc . harvard . edu /caldb / 
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the count rate in each band as described in column 6, CF is the energy conversion 
factor 0.742 cts erg'^ crri^ (and 1.837 and 0.381 for the Soft and Hard, 2-10 keV band 
respectively) appropriate for a power law spectrum with spectral index r=1.4 and 
Galactic column density Nh = 2.7 x lO^^cm^^. For sources not detected in this band, 
a 90% upper limit is reported (see Paper II for details). 

- Column 10: Full band signal to noise ratio. 

- Column 11: Full band exposure time derived from the exposure map. 

- Colurrm 12-13: The aperture photometry counts and error in the Full band (0.5- 
7 kcV) arc derived from event data for each individual Obsid and CCD where a source 
lands. Note that {Fjrate x f_exptime) 7^ fjctsjap. Circular extraction regions 
corresponding to the 90% ECF for that observation are centered on the source RA, 
Dec. The individual photometry values are then merged to produce a single set of 
values accounting for the ECF for each ObsID, given the different extraction regions 
needed. 

- Column 14: Exposure time (ksec) from the same region used to generate the aperture 
photometry. 

- Column 15-23: Same as columns 6-14 for the Soft band (0.5-2 keV). 

- Column 24-32: Same as columns 6-14 for the Hard band (2-7 keV). Fluxes and errors 
are computed for the 2-10 keV band with the conversion factor quoted above. 

- Column 33-35: Hardness ratio and 90% upper and lower errors computed as follows: 
H-S/H-l-S where H are the counts in the Hard band and S the counts in the Soft band. 
The hardness ratio was calculated starting with the EMLdetect rate values. Upper 
and lower limits were calculated using the Bayesian Estimation of Hardness Ratio 
code (BEHR, Park et al. 2006). Pseudo-source and background count values were 
generated using the net count rate, background rate (per pixel), and a 3 arcsec source 
aperture and 5-20 arcsec aperture for background areas. The aperture photometry was 
unsuitable for this purpose because the individual extraction apertures do not have the 
constant background/source area ratios required by the assumptions used in BEHR. 



5.3. Catalog Completeness &; Number Counts 

In order to provide an end-to-end check that the many calibration steps taken in deriving 
the Chandra COSMOS point source catalog have been performed correctly, we constructed 
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Table 7: Data fields in the Catalog. 



1 
1 


1\T A A/TTh' 


Chandra source name 


Z 


Source ^ 


source number. 


Q 
O 


R A 


Chandra Right Ascension (J2000, hms) 


A 

4 




Chandra Declination (J2000, dms) 


r 



pos_err 


Positional error [arcsec] 


/? 
D 


r_rate 


0.5—7 keV count rate [counts/sec] 


7 


i_rate_err 


0.5-7 keV count rate error [counts/ secj 


o 



i_nux 


0.5-10 keV thix [erg cm s J 


y 


f_flux_err 


0.5—10 keV Flux error [erg cm~^s~^] 


1 n 


f-snr 


u.o— ( Kev o/iN iiatio 


1 1 
li 


f_exptinie 


0.5—7 keV exposure time [ksec] 


1 


f_cts_ap 


0.5—7 keV aperture photometry net counts [counts] 


1 Q 


f_cts_ap_err 


0.5—7 keV aperture photometry net counts error [counts] 


1 /I 


f_exptime_ap 


0.5—7 keV exposure time from aperture photometry [ksec] 


1 c 
io 


s_rate 


0.5-2 keV count rate [counts/sec] 


lb 


s_rate_err 


0.5-z kev count rate error [counts/ sec] 


17 


sJlux 


0.5— z kev I'iux [erg cm s J 


lo 


s_flux_err 


0.5—2 keV Flux error [erg cm~^s~^] 


1 Q 


S-Snr 


u.o— z Kev o/iN xiatio 


on 


s_exptiine 


0.5—2 keV exposure time [ksec] 


91 


s_cts_ap 


0.5—2 keV aperture photometry net counts [counts] 


zz 


s_cts_ap_err 


0.5—2 keV aperture photometry net counts error [counts] 


oo 
zo 


s_exptime_ap 


0.5—2 keV exposure time from aperture photometry [ksec] 


O /I 

z4 


h_rate 


2—7 keV count rate [counts/sec] 


Of^ 
zo 


li_rat6_6rr 


2—7 keV count rate error [counts/sec] 


26 


h_flux 


2-10 keV Flux [erg cm^^s^-*^] 


27 


h_flux_err 


2-10 keV Flux error [erg cm~^s~^] 


28 


h_snr 


2-7 keV S/N Ratio 


29 


h_exptime 


2-7 keV exposure time [ksec] 


30 


h_cts_ap 


2-7 keV aperture photometry net counts [counts] 


31 


h_cts_ap_err 


2-7 kcV aperture photometry net counts error [counts] 


32 


h_exptime_ap 


2-7 keV exposure time from aperture photometry [ksec] 


33 


hr 


hardness ratio 


34 


hr_lim_lo 


hardness ratio 90% lower limit 


35 


lir_liiii_hi 


hardness rMtio 90% upper liuiil 
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Fig. 9. — The Euclidean-normalized, logN-logS curves for C-COSMOS sources with detml>W.8: 
left: Soft band (0.5-2 keV, red open circles), right: Hard band (2-10 keV). The XMM-COSMOS 
curve (black filled circles, Cappelluti et al. 2009), the soft band curve of Hasinger et al. (2005; 
green line), the Moretti et al. (2003) compilation (blue dashed line), and the CDF-N (magenta solid 
line, Alexander et al. 2003) and CDF-S (cyan solid line, Luo et al. 2008) curves. The agreement is 
good over the flux interval where the various surveys have good statistics (see text). 
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the observed logN-logS curve, i.e. the number of sources, N(>S), detected per square degree 
brighter than a a given flux, S (erg cm~^s~^) in the Soft (0.5-2 keV) and Hard (2-10 keV) 
bands. Because at the hmiting fluxes the sky coverage is smaU (Fig. 3), and so has a large 
fractional error, we used the flux limits given in Tab.l, column 3, thus omitting the faintest 
~10 sources. X-ray source counts in this flux and energy range have been well studied, giving 
us a good baseline against which to compare C-COSMOS (Cappelluti et al. 2009). 

The results, for sources detected at detml>10.8 (Tab. 5, left column) are shown in 
flgure 9, normalized by a Euclidean 1.5 slope to enable differences between various X-ray 
logN-logS curves to be seen easily. Figure 9 also shows comparisons with several other logN- 
logS curves: from Moretti et al (2003, blue dashed line), which combines data from ROSAT 
(for bright sources), XMM- Newton (for intermediate flux sources), and Chandra for faint 
sources; from Hasinger et al. (2005) logN-logS (green dashed hne); and from the CDF-N 
(magenta sohd hne, Alexander et al. 2003) and CDF-S (cyan solid line, Luo et al. 2008) 
curves. In the range where these curves overlap and C-COSMOS has good statistics the 
agreement is excellent, and C-COSMOS extends a factor ~4 below the XMM-COSMOS 
limit, as expected. 

In the Soft band, around ^2x10"^^ erg ciR-h-^, the C-COSMOS logN-logS shows a 
~20-30% underdensity at a 2(7 level with respect to the XMM-COSMOS source counts. In 
order to evaluate this deviation, we estimated the amplitude of the fluctuations expected 
due to sample and cosmic variance. According to Yang et al. (2004, 2006) and Cappelluti et 
al. (2009), the fluctuations of the counts in a box of area Q deg^ of a population of A/" deg~^ 
sources at a given flux limit, is given by a linear combination of a Poisson fluctuations and 
a cosmic variance component introduced by source clustering: 



In eq. 1 w{9) is the angular autocorrelation function expressed as a w(^) = ^ . According 
to Cappelluti et al. 2007 eq. 1 can be solved analytically by knowing the slope and the 
amplitude of w{9). By using the source surface density of Soft X-ray sources at 2x10^^"^ 
erg cm~^s~^ (i.e. ~30 source deg~^) on a box of 0.9 deg^, and assuming the angular au- 
tocorrelation function of Miyaji et al. (2007) for XMM-COSMOS (i.e. 9o^2", 7=1.8), we 
determined cr^ ~36 which corresponds to a fraction variance of 20% of the source counts. 
We can therefore conclude that a deviation of the size observed can be introduced by a 
single structure, in an area of XMM-COSMOS not covered by Chandra, that generates a 
fluctuation in the bright source counts at 1.5a level. 

Another check of the source detection efficiency at the brighter C-COSMOS flux levels 
is a comparison with the XMM COSMOS survey (Hasinger et al. 2007). As shown by 




(1) 
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Cappelluti et al. (2009) and Brusa et al. (2009, in preparation), C-COSMOS recovers ~93% 
of the XMM sources in the C-COSMOS field, resolving ~3% into close pairs. 

6. Conclusions & Future Work 

We have presented the ~0.9 sq.deg Chandra COSMOS survey (C-COSMOS) and a 
catalog of point sources from that survey. Employing a heavily overlapping tiling of ACIS-I 
observations has proven an effective method of covering a large area to a well-defined exposure 
(±12%) and uniform flux limit. The central ~0.5 sq.deg achieved an exposure of 160 ksec, 
and the outer ~0.4 sq.deg achieved an exposure of ~80 ksec. The equatorial location of 
COSMOS helped to produce a uniform tiling pattern by allowing an almost constant roll 
angle for Chandra observations over most of the target visibility window. The point source 
catalog from the C-COSMOS survey has a flux limit of 2xl0^^^erg cm^^s^^ (0.5-2 keV) and 
contains 1761 sources detected in at least one band with a probability of being spurious of 
<2xl0-^ {detml >10.8). 

The novel three-stage source detection method employed (Paper II) coped well with 
the peculiarities of the C-COSMOS tiling scheme and, more generally, is good at separating 
close pairs of sources, while retaining photometric accuracy. The C-COSMOS sky coverage 
has a sharp cut-off which produces a homogeneous flux threshold over the whole area and 
the soft band logN — logS curve for C-COSMOS matches well the Hasinger et al. (2005) 
determination over a broad flux range, giving us high confidence in the completeness of the 
catalog down to the limiting flux. 

The catalog is available in the ApJ on-line version and on the ^Chandra COSMOS 
Survey' website (see footnote 16) Supporting data products (including images, event files 
and exposure maps) are available at the ' Chandra COSMOS Survey' website and at IRSA 
(see footnote 17). 

The sub-arcsecond accuracy of the Chandra positions, together with the rich pre-existing 
deep multiwavelength coverage of the COSMOS field, allows us to reach a 96% identification 
rate for the C-COSMOS sources with counterparts in both optical and infrared, and 99.7% 
in at least one band (Paper III). 

A parallel effort on the detection of extended sources in the C-COSMOS field finds ~50 
groups and clusters (Finoguenov et al. 2009, in preparation). 

We anticipate a rich haul of science results from C-COSMOS. The Chandra sources have 
already resolved ambiguous source identifications from the XMM-COSMOS survey (Hasinger 
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et al. 2007, Brusa et al. 2007, 2008, Cappelluti et al. 2007, 2009). The paper by Fiore et 
al. (2008) on the stacking analysis of sources with extreme mid-infrared to optical ratio, 
presumably Compton Thick AGN, has been recently accepted. Several further papers are 
in preparation or submitted on: off- nuclear sources in galaxies (Mainieri et al. 2009), X- 
ray source correlation functions (Miyaji et al. 2009), the 3D cluster/AGN cross-correlation 
function (Cappelluti et al. 2009), high X-ray/optical flux ratio objects (Civano et al. 2009), 
high-redshift QSO (Civano et al. 2009), and other topics. 

A basic X-ray spectral analysis of the nearly 500 sources with more than 80 counts 
(~ 23% of the total sample) becomes possible. The resulting spectral slopes and absorbing 
column densities will allow the statistical properties of a large sample at substantial rcdshift 
and over a uniform and contiguous field to be studied effectively (Lanzuizi et al. 2009 in 
preparation) . 

There is information in C-COSMOS below the current catalog flux limit, thanks to the 
low background of Chandra ACIS. A 'stacking' analysis (Brusa et al. 2002, Hornschemeier 
et al. 2002, 2003) allows the mean X-ray properties of groups of objects to be determined. 
Miyaji et al. (2008) have solved the issues created by the C-COSMOS tihng scheme for 
stacking and papers using this tool are in preparation on z~l elliptical galaxies (Kim et 
al. 2008). The potential uses of stacking in the C-COSMOS field arc extensive, thanks to 
the multiple data sets available from which to choose samples for stacking. For example, 
there will be ~2xl0'^ galaxies with good optical spectra from z-COSMOS (Lilly et al. 2007) 
in the C-COSMOS field. This entire sample is well characterized both morphologically via 
HST imaging, and in terms of stellar population, from the UV to far-IR coverage of the 
other telescopes that have observed COSMOS (Scoville et al. 2007b). This rich data set 
will enable galaxy X-ray evolution studies by environment, morphology and luminosity using 
fine-grained stacks of C-COSMOS data with ~100 galaxies per bin, for an effective exposure 
time of ~20 Ms per bin. 

Clearly the C-COSMOS survey will be of value for some time. 
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